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ABSTRACT 

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational insta- 
bility in a collisionless plasma show that angular momentum transport by pressure anisotropy 
(P-L P\\ i where the directions are defined with respect to the local magnetic field) is comparable 
to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a 
large-scale viscosity, arises because of adiabatic invariants related to p± and p^ in a fluctuating 
magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the 
viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show 
that ~ 50 % of the gravitational potential energy is directly converted into heat at large scales by 
the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and 
magnetic energy). We show that electrons receive a significant fraction (~ [Te/Ti] 1 / 2 ) of this 
dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional 
models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow 
is greater than 0.5% for M > lO~ 4 MEdd- Thus a low accretion rate, rather than just a low 
radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. 
For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of 
3 x 10 _8 AfQyr~ 1 and an electron temperature of w 3 x 10 10 K at w 10 Schwarzschild radii; 
the latter is consistent with the brightness temperature inferred from VLBI observations. 

Subject headings: accretion, accretion disks - MHD - plasmas - Galaxy: center 



1. Introduction 

Local and global simulations have shown that 
the magnetorotational instability (MRI) is the 
likely source of turbulence and angular mo- 
mentum transpor t in many astrophysical accre- 
tion disks ( e.g.. iBalbus fc Hawlevl il99lL 1 19981: 



Hawlev. Gammie. fe Balbus ll995HArmitage Ill998f) 



The MRI has been studied in a wide variety 
of contexts, from t h e disk of the Milky Way 



( Piontek fc Ostriker 



20051) to planet-forming 
disks around young stars (e.g., Winters, Bal- 
bus, & Hawley 2003). Our recent work has fo- 
cused on the e volution of the MRI in collision- 
less plasmas (jQuataert. Dorland. fc Hammett 
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2002, Sharma Ha mmett, & Quataert 2003, 
Sharma et al. 20061 [here after SHQS]; for related 



work, see Balbus 1 120041 and Krolik & Zweibel 
20061) . These calculations are motivated by the 



application to radi atively inefficient accretion flow 
(RIAF) models (IchimaruJ 119771 ; iRees et. al. 



19821 : iNaravan fc Yi I Il995l ). which are believed 



to be applicable to accretion onto black holes and 
neutron stars when the accretion rate is less than 
a few % of Eddington. Sgr A* in our Galactic 
Center is the paradigm ca se for accretion via a 



RIAF (see iQuataert 1120031 for a review). 



One of the central problems in the theory of 
hot accretion flows remains understanding the ex- 
treme low luminosity of accr etion onto black holes 



in m any galactic nuclei (IFabian fc Canizares 



19881 IFabian fc Reesl I1995U Loewenstein et al 



20011 ) and X-ray binaries (e.g., Narayan, Mc- 
Clintock, & Yi 1996). In the context of RIAF 
models, the observed low luminosities can be ex- 
plained if most of the available mass is not ac- 
creted because of outflows or convection (Bland- 
ford & Begelman 1999; Quataert & Gruzinov 
2000; Stone, Pringle, & Begelman 1999; Narayan, 
Igumenshchev, & Abramowicz 2000), or if the 
gravitational potential ener gy primarily heats 



the poorly radiatin g ions (jRees et. al. 1982; 



Naravan fc Yilll995f ). Global MHD simulations 



support the former possibility Jstone fc Pringle 
120011 : lHawlev. Balbus. fc Stone Il200ll ). 

A number of analytic arguments have been pre- 
sented for the heating of electrons and ions in the 
collisionless plasmas appropriate to RIAFs. These 
estimates focus on heating by MHD turbulence 



(IQuataert1l998t[Gruz 



IQuataert fc Gruzinov 



nov1ll998l:lBlackman 11199 
19991: iMedvedevI bo~ 



or r econnection ( Bisnovatvi-Kog an fc Lovelace 



19971 : IQuataert fc Gruzinov 1 1 19991 ). Our shearing 



box simulations of MRI turbulence in a collision- 
less plasma suggest an additional mechanism for 
particle heating in RIAFs. Adiabatic invariance 



(/; 



/B oc p±/B = constant) in a collisionless 



plasma with a fluctuating magnetic field results in 
a macroscopic pressure anisotropy, with the pres- 
sure perpendicular to the field lines {p±) exceed- 
ing that along the field lines (p||). This pressure 
anisotropy results in a viscous transport of angu- 
lar momentum that is comparable to the trans- 
port by magnetic stresses (SHQS). In a collisional 
plasma, the magnitude of pressure anisotropy, and 



thus the viscosity, is determined by Coulomb col- 
lisions. In turn, the relative heating of electrons 
and ions by viscosity is set by the temperature and 
particle-mass dependence of the Coulomb cross- 
section (which results in primarily ion heating). 
In a collisionless plasma, however, the pressure 
anisotropy — and thus the viscosity — is regulated 
by the growth of small scale instabilities that vio- 
late adiabatic invariance and isotropize the plasma 
pressure (e.g., the firehose, mirror, ion cyclotron, 
and electron whistler instabilities). As we shall 
show, these instabilities also regulate the viscous 
heating of electrons and ions in a collisionless ac- 
cretion flow. 

In §2 we summarize the equations that describe 
the evolution of the pressure tensor (pn and p±); 
a more detailed presentation is given in SHQS. 
We then present upper limits on the pressure 
anisotropy because of gyroradius scale instabili- 
ties. We use these limits to provide analytic es- 
timates for the viscous heating of electrons and 
ions in a collisionless plasma and contrast these 
results with the more familiar results appropriate 
to a collisional plasma. In §3 we discuss the results 
of two-component (electron + proton) shearing- 
box simulations of the MRI. This is an extension 
of our single component simulations (SHQS). In 
§4 we use 1-D models of RIAFs to calculate the 
radiative efficiency using the viscous heating of 
electrons calculated in §2 and §3. We also show 
that electron Coulomb collisions are less effective 
at reducing the pressure anisotropy in RIAFs than 
kinetic microinstabilities and thus that the kinetic 
approach focused on throughout this paper is ap- 
propriate. Finally, §5 is a discussion of our key 
results and their application to RIAFs. 

2. Analytic considerations 

In the limit that the fluctuations of interest 
have length scales much larger than the gyrora- 
dius and time scales much longer than the cy- 
clotron period, the equations des cribing the evo- 
lution of jP || and p 1 are given by dKulsrud 



Snvder. Hammett. fc Dorland 



19971 



and refer- 



ences therein) 



dp\\, s 
dt 



+ V-(p||, a V)+V-q||, a + 2p|| )S bb: W 



2q ±!S V • b = -^eff,s(Pls -P±,s)(l) 



2 



dp±,s 
dt 



+ V -(px,.V) + V -q±, s +p±, s V -V 
- p±, s bb : VV + q ±!;5 V • b 

= -q^Z/.sCpi-.a (2) 



where the subscript 's' denotes the species (elec- 
trons and ions [we only consider protons] ) , b is the 
local magnetic field direction, V is the fluid veloc- 
ity, = <?j_b (q|| = giib) is the flux of perpen- 
dicular (parallel) thermal energy along the local 
magnetic field. The effective pitch-angle scatter- 
ing rate due to high frequency waves (or Coulomb 
collisions) is given by v e f / . 

Equations (JTJ) and © can be obtained by tak- 
ing the moments of the drift kinetic equation, the 
Vlasov equation in the limit of large wavelength 
and long timcscales. These moment equations run 
into the usual "closure" problem. In our numerical 
calculations, we use a simple closure for the heat 
flux, in which the heat flux is proportional to the 
temperature gradient, with the conductivity being 
a free parameter. This model can roughly repro- 
duce kinetic effec ts such as collisionless damping of 
linea r modes (s ee Sharma. Hammett. fc Quataert 



2003 : Sharma 2006^ 1 . For the analytic estimates 
in this section, these details of the heat fluxes are 
not important. In addition, our numerical simu- 
lations show that the electron heating rate does 
not depend strongly on the details of the heat flux 
model. 

Equations (J) & © can be combined to derive 
equations for the internal energy (e s = 3p s /2 = 
pii s /2 + p±. s ) and pressure anisotropy (Ap s = 
P\\. s - P±,s), which are given by 



de s 
dt 

dA Ps 
dt 



V- (e s V + q s )+p_L, s V- V 
-Ap s bb : VV, (3) 
V- (A Ps V + Aq s ) = -3p s bb : VV 



+ P±,sV ■ V + 3q ± , s V ■ h 
- v ef f iS Ap s , 



(4) 



where q s = q||, s /2 + q_L, s and Aq s = q^ a - q ±s . 
In statistically steady MRI-driven turbulence, 
Ap s < on average (see, e.g., Fig. 4 of SHQS). 
This sign of the pressure anisotropy corresponds 
to the expected outward transport of angular mo- 
mentum by the viscous stress. 



If the fluctuations of interest are relatively in- 
compressible then the volume averaged version of 
equation ([3]) is given by 



3_d 

2dt Ps 



4, 



-Ap s bb : VV. 



(5) 



The RHS of equation ([5]) represents heating due 
to the anisotropic pressure. In an accretion 
flow, the velocity gradient in equation ([5]) can be 
roughly decomposed into two components, that 
due to the background shear O(r), and that due 
to the turbulent velocity fluctuations SV (such 
a decomposition is formally made in the shcar- 
ing box simulations described in §3, see e.g., 
Hawlev. Gammie. fc Balbus I Il995l ). Thus the 
heating rate itself can be decomposed into two 
contributions: a heating term given by 



din 7' 



Ap s b r b,f 



and one given by 

q+ 2 s = -Ap s bb : VSV. 



(6) 



(7) 



Equation ([6]) represents viscous heating due to the 
background shear while equation ([7]) represents 
dissipation of the turbulent fluctuations (by, e.g., 
collisionless Landau damping) 

In the analytic estimates that follow we focus 
on viscous heating due to the background shear; in 
our simulations, we find that this is the dominant 
contribution to the total heating in equation l[5]). 
Note that, just as for a collisional viscosity, this 
term represents the direct conversion of ordered 
kinetic energy (differential rotation) into heat at 
large scales o The other sources of heating are 
less direct: energy is converted into magnetic and 
velocity fluctuations which are dissipated by colli- 
sionless damping at large scales, or at small scales 
via a nonlinear cascade. 

We now consider the collisional (Braginskii) 
and collisionless limits of equation (O and quan- 
tify the viscous heating in each of these limits. 



our numerical simulations, there is also energy lost due to 
grid-scale viscosity and resistivity; these will be discussed 
in §3. 

2 Collisionless viscous heating is analogous to the betatron 
mechanism (also known as magnetic pumping) where adia- 
batic invariance and pressu re isotropizing collisions can re - 
sult in net plasma heating ( Budkcr 1961; Kulsrud 2005). 
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2.1. The Collisional limit 



2.2.1. Upper Limits on the Pressure Anisotropy 



If the collisional mean free path is small com- 
pared to the gradient length scales of interest 
then the viscous stress in a plasma can be writ- 
ten as a pressure tens or, anisotropic w ith re- 



spect to the field lines (|Braginskii I Il965h . where 



the magnitude of the pressure anisotropy is lim- 
ited by collisional pitch angle scatterin g (see, 
e.g., ISnvder. Hammett. fe Dorland I Il997l for the 
equivalence of Braginskii viscosity and anisotropic 
pressure). Pressure anisotropy in the colli- 
sional limit of equation ((U is given by Ap s = 
— (p s /iA;)3bb : VV, where v s is the (Coulomb) 
collisional pitch angle scattering rate. Using 
3 ^ 2 equation ((5]) becomes 



-1/2™ 



1v, 



oc m 



V 2 T^ 2 (bb 



VVJ 



(8) 



where m s and T s are the mass and temperature 
of the species, respectively. As noted above, for 
an accretion flow bb : VV ~ dil/d\nr. Equation 
([8]) implies that for T e = Ti, the heavier ions are 
heated preferentially by viscosity by a factor of 
~ 40. This conclusion is strengthened if T e < Ti, 
as is typically the case in hot accretion flows. 

2.2. The collisionless limit 

When the mean free path is comparable to or 
larger than the gradient length scales of interest, 
the collisional (Braginskii) description is no longer 
valid. In the collisionless regime, the pressure 
anisotropy is limited by pitch angle scattering due 
to cyclotron frequency fluctuations that violate 
adiabatic invariance. A minimum level of pitch 
angle scattering is that due to kinetic instabilities 
generated by the pressure anisotropy itself. We 
focus on that possibility here. In principle, the 
rate of pitch angle scattering can be much higher 
if high frequency fluctuations are efficiently gen- 
erated by a turbulent cascade, shocks or recon- 
nection. Whether this in fact occurs is difficult 
to assess and may depend on the specific prob- 
lem of interest; theories of incompressible MHD 
turbulence, however, imply little power near the 
cyclo t ron frequency (e.g ., iGoldreich fc Sridhar 



1995t iHowes et al.tl2006li . and measurements h 



the s olar wind show significant pressur e anisotropy 
(e.g., Kasper. Lazarus, fe Gary 2002h . 



To quantify the rate of pitch angle scattering 
in a collisionless plasma, we first discuss limits 
on the pressure anisotropy due to kinetic insta- 
bilities. As discussed in §1 (and SHQS), pressure 
anisotropy is naturally created in MRI turbulent 
plasmas due to the fluctuating magnetic field. For 
p± > pn, as occurs if the magnetic field is am- 
plified by the MRI, the relevant instabilities are 
the mirror and ion-cyclotron instabilitie s for the 
ions (jrlasegawa Hl969UGarv fc Lee 1994; also dis- 
cussed in SHQS), and t he electron whistler in sta- 
bility for the electrons (jGarv fc Wang 1 11996) . In 
the opposite limit (pu > p\_), the fire hose instabil 



ity reduces the pressure anisotropy (jGarv et al. 
19981 : iLi fc Habbal Il2000l SHQS). 



The mirror instability arises when the pressure 
anisotropy becomes larger than p±,i/p\\,i — 1 > 
l//3_L,j. However, the fastest growing mode occurs 
at the gyroradius scale (and hence violates /i in- 
variance and leads to pitch angle scatt ering) only 
when (SHQS; see also iHasegawa 1119691 ) 



Phi 



1 > 



7 



(9) 



For a pressure anisotropy smaller than this thresh- 
old, no significant pitch angle s cattering occurs 



(jMcKean. Gary, fc Winske~lll993f) 



A plasma is formally unstable to the ion cy- 
clotron instability for arbitrarily small values of 
T± i i/Tu j— 1 > 0. However, the instability can only 
grow on a reasonable timescale (~ rotation period) 
when the pressure anisotropy exceeds a threshold 
which can be written as (e.g., iGarv fc Lee 19941 
SHQS) 

^ (10) 

i ll,< P\\,i 

where Si and oti depend on 7/^, the ratio of the 
growth rate to the ion cyclotron frequency. We 
have used the pub lically available WHAMP code 
(jRonnmark 1982) to calculate Si and «i for a 
specified value of j/fli. For RIAFs, fli ^> O, 
the disk rotation frequency. Because pressure 
anisotropy is generated on the turnover time of 
the turbulent fluctuations, which is comparable to 
the rotation period of the disk (~ S! -1 ), 7 
is sufficient to generate significant high frequency 
fluctuations which provide pitch angle scattering. 
Taking = 10 -4 as a fiducial estimate of 
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the growth rate required for significant pitch an- 
gle scattering, we find Si — 0.35 and on = 0.45 
as the threshold for the ion cyclotron instability. 
For 7 <C f2j the pressure anisotropy threshold (Si) 
depends very weakly (roughly logarithmically) on 
the growth rate so our results are not sen- 

sitive to the parti cular choice of j/Cli (see, e.g., 
Gary fc Leelll994h . 

As noted above, the mirror instability and ion 
cyclotron instability both act to isotropize the 
plasma when p± > p^ . For accretion disks with 
0i ^ 100, the ion cyclotron instability is gener- 
ally more important than the mirror instability. 
In our analytic estimates we assume that the ion 
cyclotron instability dominates, but our simula- 
tions include both the mirror and ion cyclotron 
thresholds discussed above. 

Electrons with an anisotropic pressure (p_i_, e > 
P\\, e ) are unstable to the whistler instability, which 
has a frequency ~ VL e and thus violates the adia- 
batic invariance of the electrons, producing pitch 
angle scattering. The threshold for the growth of 
the whistler instability with j/Cl e = 5 x 10~ 8 is 



given 



1 > 



Se 



(11) 



with ,!5 e = 0.13 and a = 0.55 (again using the 
WHAMP code) . As with the ion cyclotron thresh- 
old, the whistler threshold depends only weakly on 
the choice of j/Fl e - 

In our simulations, we find that p±. s > p\\ lS 
holds over most of the computational domain, as 
expected because of the outward transport of an- 
gular momentum by the viscous stress. There are, 
however, regions where the pressure anisotropy 
has the opposite sign. In this case electron and 
ion firehose ins tabilities wi l l limi t the pressure 
aniso tropy to ([Gary et al. I Il998l iLi fc Habbaf 
2000h 



1 



< 



(12) 



The magnitude of the viscous heating in equa- 
tion ([5]) is directly proportional to the pressure 
anisotropy Ap s , which is < on average and sat- 



3 For non-relativistic electrons, this choice for the electron 
growth rate corresponds to the same value of 7 as that 
used in determining the ion cyclotron threshold for the ions 
above (f/Sli = 10~ 4 ). 



urates at a value close to the p± tS > pu s thresh- 
olds discussed above (as shown explicitly in Fig. [3] 
discussed in the next section). For an analytic es- 
timate, we assume that the ion and electron pres- 
sure anisotropics are bounded by the ion cyclotron 
and whistler thresholds, respectively. We then es- 
timate the viscous heating of electrons and ions 
in the collisionless limit using equations ([5|), (|10D . 
and fTTj) 

g+.ocS.^ 1 - ". (13) 

The ion cyclotron and whistler thresholds corre- 
spond to Si/5 e w 3 and oti ~ a e ss 0.5, in which 
case qy tx T^ 2 . Viscous heating in a collisionless 
plasma is thus quite different from that in a colli- 
sional plasma: the electron and ion heating rates 
are comparable when T e ~ Ti. Indeed, in the ab- 
sence of radiative cooling, equations (| 13|) and (J5]) 
imply that the electron to ion temperature ratio 
approaches T e /Ti ~ (S e /Si) 2 ~ 0.1 at late times, 
if viscous heating is the dominant heating mech- 
anism (even if the ions are initially much hotter 
than the electrons). 

3. Numerical Simulations 

In our previous work (SHQS) we modifie d the 



ZEUS MHD code (jStone fc Norman I Il992albl ) to 



model the dynamics of the MRI in a collisionless 
plasma, by including pressure anisotropy, ther- 
mal conduction along field lines, and subgrid mod- 
els for pitch angle scattering by microinstabilities 
when the pressure anisotropy exceeds the thresh- 
olds (discussed in §2.2.1). Our previous work mod- 
eled a single fluid, effectively the ions since they 
typically dominate the pressure if 3> T e . We 
now extend this work to include both ions and 
electrons. As in SHQS, we restrict ourselves to 
local shearing box simulations. 

We solve the mass and momentum conservation 
equations for the shearing box in MHD (eqs. [35] 
and [36] in SHQS), where the pressure is the sum 
of electron and ion pressures, each evolved accord- 
ing to equations (P) and ^ . Heat fluxes parallel to 
the field lines (q_i_ and qii ) are given by equations 
(40)-(44) of SHQS applied to electrons and ions. 
The conductivity along the field lines is given ap- 



proximately by k 



'in 



/{vthkL 



where v t h 



is the thermal velocity of the species, 24/7 is the 
pitch angle scattering rate due the microinstabil- 
ities, and &l is a parameter that corresponds to 
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the typical wavenumber for Landau damping; this 
gives the correct collisionless damping rate for a 
fluctuation whose wavenumb er is kh but is only 
approximate for other modes ( Sharma 2006). For 
small feffj the conductivity in the heat flux is 
oc fc^ 1 so that small values of Ul correspond to 
rapid heat conduction while — > oo is the adia- 
batic (CGL) limit with no heat conduction^] 

In our simulations, the usual shearing box 
boundary conditions are used (Hawley, Gammie, 
& Balbus 1995). Initial conditions are described 
in SHQS. The present simulations use a purely 
vertical initial field with ft = 400. In all runs we 
initiate a specified electron and ion temperature 
ratio at the end of 5 orbits so that our results 
are not affected by the strong channel phase at 
~ 4 orbits that represents the initial non-linear 
saturation of the MM. Electron and ion cooling 
is neglected. Unless specified otherwise, the res- 
olution is N x x N y x N z = 27 x 59 x 27, and 
the box size L x x L y x L z = 1 X 2n X 1. The 
initial parameters in code units are: ion pressure 
Po = 10~ 6 , rotation rate £1 = 10~ 3 , and initial 
density p = 1 (as in SHQS). The thermal conduc- 
tion parameter = 0.5/Sz for both electrons and 
ions, where Sz (= L z /N z ) is the grid spacing in 
the vertical direction. Note that this value of the 
conductivity corresponds to correctly capturing 
Landau damping for modes with wavelengths of 
(4ir/N z )L z w 0-5L z , i.e., for modes whose wave- 
lengths are approximately half the size of the box. 

Since ZEUS is not a conservative code, energy is 
lost because of numerical dissipation. However, to 
quantify the importance of different heating mech- 
anisms it is desirable to quantitatively track the 
energetics of the accretion fl ow. To do so, we have 

" ( 20031) 



implemented the method of iTurner et al. 



in which energy conservation in ZEUS is signifi- 
cantly improved by keeping track of the kinetic 
and magnetic energies lost because of grid-scale 
averaging. This energy can also, if desired, be 
added as a source of plasma heating. However, we 
do not do so because it is physically unclear which 
species should receive this energy. In what follows 
we identify energy lost numerically in the magnetic 
field update as "numerical resistive" heating (g^ft) 
and energy lost in the transport step for momen- 



4 CGL refers to the double adiabatic model of 
IChew. Glodberger. fc Lewi dl956h . 



turn as "numerical viscous" (q^ v ) heating. Com- 
pressive heating, — p±\7 ■ V, is negligible because 
the MRI is relatively incompressible. In addition 
to numerical heating due to grid-scale dissipation, 
there is also viscous heating at large-scales due 
to pressure anisotropy {qy\ see eq. [5]). This is 
captured by our simulations — including the cor- 
rect ratio of electron to ion heating — and repre- 
sents direct conversion of gravitational potential 
energy into heat at large scales (see §2). Figure 
Q] shows a plot comparing the change in the total 
plasma energy with the work done by the bound- 
aries in a typical shearing box simulation. Energy 
is conserved to within 3 x 10 -3 over a disk rotation 
period, which is adequate for our purposes. 

3.1. Stresses and heating rates 

Tables [T] & [5] summarize the properties of a 
number of simulations with the ion to electron 
temperature ratio initialized at values from 1 — 10 4 
after the MRI saturates at 5 orbits. We find 
that the properties of saturated MRI turbulence 
in a collisionless plasma are relatively insensi- 
tive (to within a factor of few) to the conduc- 
tivity (parameterized by Ul) or numerical resolu- 
tion. In particular, angular momentum transport 
quantities like the anisotropic stress {a a, where 
a is the stress divided by the initial pressure), 
Maxwell stress (q:m), and Reynolds stress {olr), 
and the associated viscous (qy), "numerical resis- 
tive" (q^fj)-, and "numerical viscous" {q% v ) heat- 
ing rates do not depend sensitively on k^ or res- 
olution. Instead, the most important physical ef- 
fect in the evolution of the MRI in a collisionless 
plasma is pitch angle scattering by kinetic insta- 
bilities (§2.2.1), which determines the magnitude 
of the anisotropic stress and thus the magnitude 
of the viscous heating. For the pitch angle scat- 
tering thresholds described in §2.2.1 we find that 
anisotropic stress and the associated viscous heat- 
ing are the dominant terms in the angular momen- 
tum and internal energy equations. In addition, 
viscous heating is dominated by direct heating via 
the large-scale shear (eq. [BJ; qy^; heating due to 
the turbulent fluctuations (eq. [7J; qy 2 ) is signifi- 
cantly smaller (see Tables [1] & [2]) . 

Figure [2] shows the volume average of different 
heating rates as a function of time in our shearing 
box simulations with Ti/T e = 10 4 at 5 orbits, for 
N x x N y X N z = 27 x 59 X 27 and k L = 0.5/Sz 
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(left), and N x x N y x N z = 54 x 118 x 54 and 
k,L — 0. 25/(5 z (right). These combinations of reso- 
lution and kh correspond to the same conductiv- 
ity in the two different simulations. Both the high 
and low resolution simulations show statistically 
similar heating results. Direct viscous heating is 
slightly larger than numerical resistivity and vis- 
cosity, implying that ~ 50% of the gravitational 
potential energy of the accretion flow is directly 
converted into heat at large scales (see also Table 
[1]) . Figure [2] also shows that the different heat- 
ing rates are correlated and fluctuate together in 
time. Although the viscous heating rate varies as 
a function of time, we find that the ratio of the ion 
to the electron viscous heating does not show as 
large statistical variations, i.e., the ion and elec- 
tron heating rates increase/decrease together in 
tandemd 

Figure [3] shows the volume averaged pressure 
anisotropy for electrons and ions with Ti/T e set 
to 10 at the end of 5 orbits; the mirror, ion cy- 
clotron, and electron whistler thresholds are also 
shown. The ion pressure anisotropy is larger than 
that of the electrons for two reasons: first, the 
electron pressure anisotropy threshold is smaller 
than that of the ions (S e ~ OASi), and second, 
the electron (3 is smaller by a factor of 10. Figure 
[3] also shows that the electron pressure anisotropy 
is relatively close to the electron whistler instabil- 
ity threshold and that the ion pressure anisotropy 
roughly tracks the ion cyclotron threshold. 

In a collisionless plasma, the MRI grows faster 
for an initial magnetic field configuration with 
Bcj) = B z than for = (Quataert, Dorland, 
& Hammett 2002; Sharma, Hammett, & Quataert 
2003). However, the nonlinear saturated state 
is similar to th at of the pure vertical field case 



described here ( Sharma 20061 ). In particular 



the anisotropic stress is again comparable to the 
Maxwell stress and viscous heating accounts for 
a significant fraction of the total heating of the 
plasma. 

To quantify the relative heating of ions and elec- 
trons as a function of Ti/T e , we initialize simula- 
tions in the saturated turbulent state (after 5 or- 



5 This is because the ion and electron pressure anisotropics 
are approximately equal to the ion cyclotron (for ions) and 
electron whistler (for electrons) thresholds, which have a 
similar dependence on f3 (see eq. |13j). 



bits) with T,/T e = 10 4 ,10 3 ,10 2 ,10, & 1. Figure 
|4] shows the ratio of the volume averaged ion and 
electron temperatures for kh = 0.5/ Sz, and in the 
CGL limit (q± = = 0). Initially cold electrons 
become heated significantly by viscous heating as 
discussed in §2.2. Based on extrapolating the re- 
sult of Figure [4] to even later times, we estimate 
that the late-time value of Ti/T e due to viscous 
heating alone is ~ 10 — 30, in reasonable agree- 
ment with the analytic estimate in §2.2. 

Figure [5] shows the ratio of the volume aver- 
aged ion and electron heating rates as a function 
of Tj /T e for the same calculations as in Figure |4] 
('o' for CGL, '+' for k L = 0.5/ 'Sz, and 'A' for a 
higher resolution simulations with k^ = 0.25/ Sz). 
This plot was made by averaging the heating and 
Ti/T e over 0.1 orbits in a number of different simu- 
lations (so that the temperature ratio is fairly con- 
stant over the time of averaging) . The middle solid 
line in Figure [5] shows q^Jq^ e = 3(T 4 /T e ) 1 /2 j 
which is roughly the analytic prediction from §2.2 
(eq. [13] ) assuming that the ion and electron pres- 
sure anisotropics are set by the ion cyclotron and 
whistler instabilities, respectively. The agreement 
between the analytic prediction and the numeri- 
cal results is particularly good in the absence of 
conduction: the heating ratio Qyi/lve 1S slightly 
larger (less than a factor of two) with conduction 
than without it. Figure \5\ also shows that the elec- 
tron to ion heating ratio does not depend sensi- 
tively on the resolution of the simulation. 

4. Implications 

To quantify the importance of viscous heating 
for observations of accreting black holes, we use 
a heating prescription motivated by our analytic 
and numerical resul ts in one-dimensional model s 
of RIAFs (based on Quataert fc Naravan 1999j) . 
From equation (|13p and Figure [5j we approxi- 
mate the fraction of the viscous energy that heats 
the electrons as * e q+/q+ w 0.33(T e /T l ) 1/2 . 
This approximation is consistent with our numer- 
ical simulations in the CGL limit; 5 is slightly 
smaller if thermal conduction is included (Fig. 
[5]). Our one-dimensional calculations include elec- 
tron cooling by free-free emission, synchrotron 
radiation, and Inverse Compton emission. The 
electron temperature, the spectrum of radiation, 
and the radiative efficiency are calculated self- 
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consistently given the assumed electron heating 
rate and the density, radial velocity, etc. from a 
one-dime nsional dynamical model. Fo r more de- 
tails, see iQuataert k, Naravan I (jl999h and refer- 
ences therein. 

Figure [6] shows the resulting radiative efficiency 
of R1AF models as a function of accretion rate; we 
assumed a = 0.1 and /? = 10 in these calculations 
but the results are not that sensitive to reason- 
able variations in these parameters. Three choices 
of 5 (S = OmlTjTi] 1 / 2 , S = 0.33PVT,] 1 / 2 , and 
S = O.nlTe/Ti] 1 / 2 ) are shown, given the uncer- 
tainty in the exact magnitude of the viscous heat- 
ing (e.g., the relative contribution of viscous heat- 
ing and other heating mechanisms, the effects of 
thermal conduction, and the uncertainty in the 
pressure anisotropy thresholds). For concreteness, 
the calculations shown in Figure [6] assume that the 
accretion rate is constant with radius. Since most 
of the radiation is produced at small radii, the 
radiative efficiency only depends on the accretion 
rate at ~ 10 Schwarzschild radii and so the x-axis 
in Figure [6] can be interpreted as such. Figure [6] 
demonstrates that with analytically and numer- 
ically motivated viscous heating rates, the radia- 
tive efficiency is > 0.5% for M > 10~ 4 M Edd for all 
of our electron heating models considered. At very 
low M, the efficiency decreases in all models be- 
cause the electron cooling time (primarily by syn- 
chrotron radiation) becomes longer than the inflow 
time in the accretion flow. Since our calculations 
in Figure [6] do not account for the ~ 50% of the 
gravitational potential energy lost via grid-scale 
averaging in our shearing box simulations (see §3), 
the resulting radiative efficiencies are likely a lower 
limit to the true radiative efficiency. 

4.1. Electron collisionality 

One possible caveat to the application of our 
fully collisionless results to RIAF models is that 
if the electron Coulomb collision frequency v e is 
significantly larger than the rate at which pres- 
sure anisotropy is created by turbulence, then 
Coulomb collisions will suppress the electron pres- 
sure anisotropy faster than kinetic instabilities^ 
As a result, electron heating by anisotropic pres- 



3 Ion isotropization is much slower than electron isotropiza- 
tion and so we do not need to consider a similar effect for 
the ions. 



sure will be negligible. We thus estimate an up- 
per limit on the mass accretion rate above which 
collisions will be able to wipe out the pressure 
anisotropy. This calculation is analogous to stan- 
dard estimates in the literature of the critical 
mass accretion rate ~ a AfEdd above which it 
is not possible to maintain a two-temperature 
RIAF because Coulomb collisions thermally cou- 
ple the electrons a nd protons on an inflow time 
(jRees et. al. II1982I) . Because electron pitch angle 
scattering is ~ m p /m e ~ 10 3 faster than electron- 
proton energy exchange, the critical accretion rate 
above which the electrons are isotropized might 
be expected to be ~ 10~ 3 a 2 AfEdd- This result is 
not correct, however, for two reasons: 1. Pressure 
anisotropy is created on a timescale even shorter 
than an eddy turnover time due to the fluctuat- 
ing magnetic field, rather than on the inflow time 
(see eq. [15] discussed below). 2. The electron 
Coulomb collision frequency depends strongly on 
electron temperature, which is itself a strong func- 
tion of accretion rate. 

The average pressure anisotropy can be esti- 
mated from equation ([4J. The second term van- 
ishes after spatial averaging. Assuming that fluc- 
tuations are incompressible, and neglecting the 
heat fluxes, the dominant terms are 



d_ 

dt 



3p e bb : VV -v e f f, e Ap e -v e Ap e (14) 



where we have added a separate term for the 
rate v e of pitch angle scattering by Coulomb col- 
lisions to clearly distinguish between the effects of 
Coulomb collisions and microinstabilities (^ e //, e ). 
If v e is negligible, then steady state occurs when 
the pitch angle scattering caused by microinstabil- 
ities, at the rate f e //,e balances the first term on 
the RHS of equation (| 14[) . The scattering by mi- 
croinstabilities occurs very rapidly (on the electron 
cyclotron timescales) if the pressure anisotropy is 
larger than the threshold to excite the instabili- 
ties, so in practice v B ff,e will always be just large 
enough to keep Ap e near the threshold. Thus bal- 
ancing the first and second terms on the RHS gives 



3 J^bb : VV 



3517/3 



1/2 



1/2 



(15) 

using a e and S e from §2.2.1. If the collisional 
pitch angle scattering rate exceeds that estimated 



S 



in equation (|T5|) , then electron viscous heating will 
be reduced relative to the pure colhsionless calcu- 
lations presented in §2 & 3. Figure [7] shows this 
explicitly by plotting the average electron heat- 
ing rate for shearing-box simulations that include 
an additional Coulomb collision rate in addi- 
tion to the subgrid models for kinetic instabilities. 
The numerical results in Figure [7] show that above 
a critical v e , electron viscous heating is substan- 
tially reduced; this critical v e depends on T e /Ti, 
in good agreement with the analytic estimate in 
equation (fT5|) . 

To apply the above results to RIAF models, we 
note that for the accretion rates of interest, the 
electrons are only marginally relativistic (see Fig. 
[8] discussed below) and so we use the Coulomb col- 
lision frequency for non-rel ativistic electro n s. Us- 
ing the ADAF solution of lNaravan fc Yi I (|l995h 



to estimate the density and ion temperature in the 
accretion flow, we find that electron Coulomb col- 
lisions are negligible for the isotropization of the 
electron pressure for accretion rates satisfying 

where r is the distance from the central object 
in units of the Schwarzschild radius and 9 e = 
kT e /m e c 2 is the electron temperature normalized 
to the electron rest mass energy. 

Figure [5] shows the electron temperature as a 
function of radius for our one-dimensional RIAF 
models with 5 = O^Te/T;) 1 / 2 , for a variety of 
accretion rates (solid lines). The protons are at 
roughly the virial temperature for all M since they 
do not cool significantly. By contrast, the elec- 
tron temperature depends strongly on M, with 
the temperature decreasing at small radii for in- 
creasing M. For M < 10 _5 MEdd, however, the 
electron temperature is relatively independent of 
M because the electron cooling time becomes long 
compared to the inflow time. The dotted line 
shows the electron temperature for low M for 
S = 0.66(T e /T i ) 1 / 2 ; the temperature is a factor 
of ~ 2 larger at small radii in this case because of 
the additional heating. 

Figure [5] shows that 9 e ~ 1 at M ~ lCr 2 M Edd ; 
at this high M, the electron temperature is rela- 
tively independent of our assumptions about the 
electron heating. This is because the electrons are 
efficiently cooled by synchrotron radiation and In- 



verse Compton emission, and become marginally 
relativistic. Using 8 e ~ 1 in equation (|16p. we con- 
clude that electron Coulomb collisions are unim- 
portant for isotropizing the electrons for M <C 
0.06MEdd- Since the electron temperature in- 
creases with decreasing M (Fig. [8]) , Coulomb colli- 
sions rapidly become unimportant for isotropizing 
the electrons below ~ 0.06A/Edd- We thus con- 
clude that the colhsionless calculations of the ra- 
diative efficiency in Figure [5] should be applicable 
for essentially all of the accretion rates considered. 

5. Discussion 

We have shown that viscosity arising from 
anisotropic pressure is a significant source of heat- 
ing in hot accretion flows (RIAFs). In shearing 
box simulations of the MRI in colhsionless plas- 
mas, viscous heating is comparable to the rate 
at which energy is lost by grid-scale numerical 
averaging. It thus accounts for ~ 50% of the 
gravitational potential energy released by the in- 
flowing matter. In a real accretion disk, the fate 
of the remaining ~ 50% of the energy that is dissi- 
pated at small scales is unclear. Although viscos- 
ity along the field lines due to pressure anisotropy 
can damp parallel motions, it cannot dissipate en- 
ergy in motions perpendicular to the field lines 
(e.g., Alfven modes). This energy is presumably 
dissipated through colhsionless damping of fluctu- 
ations in a turbulent cascade at the ion L armor 
radius (jQuataert 111999 ; iHowes et al. I l2006h . 

In contrast to the kinetic and magnetic en- 
ergy lost to grid-scale averaging, the physics of 
viscous heating at large scales is well-captured 
by MRI simulations and can be accurately ap- 
proximated using a simple analytic expression 
that depends primarily on the average pressure 
anisotropy in the plasma (eq. |13|V In turn, 
the magnitude of the pressure anisotropy is set 
by pitch angle scattering due to small-scale ki- 
netic instabilities (see §2.2.1). We thus conclude 
that pressure anisotropy plays an essential role in 
both the energy (via viscous heating) and angu- 
lar momentum (via anisotropic stress) budgets of 
RIAFs. This interplay of pressure anisotropy and 
microinstabilitics is likely to be important in many 
other weakly collisional plas mas, e.g., the solar 



wind and X-ray clusters (e.g.. lSchekochihin et al 
20051 ). The only way out of this conclusion is 
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if there is efficient generation of high frequency 
fluctuations that provide pitch angle scattering 
(and thus pressure isotropization) at a rate much 
faster than that of the kinetic instabilities dis- 
cussed in §2.2.1. This does not appear to be true 
in the solar wind, whe re large pressure anisot ropics 



are measured (e.g., Cranmer et al. 1999f ) and 



where there is in situ experimental evidence for 
pressure isotropization via the fi rehose instability 
(jKasper. Lazarus, fc Gary 1120021 ). 

In addition to viscous heating, numerical simu- 
lations can in principle also capture heating by the 
collisionless (Landau/Barnes) damping of large- 
scale turbulent fluctuations. In particular, the 
slow magnetosonic mode is strongly damped in 
a collisionless plasma even on large scales and 
has been proposed as a source of proton heat- 
ing m RIAFs (e.g., iQuataert 1 1199a iBlackman I 
1999) We find little direct evidence for such 



heating in our numerical simulations; the part of 
the anisotropic pressure heating (eq. [5]) propor- 
tional to the background shear dfl/dlnr (eq. [6]) 
dominates over the part proportional to the tur- 
bulent velocity fluctuations (eq. [7]) for all of our 
simulations (see Table [T]). In addition, the rela- 
tively weak dependence of the plasma heating on 
the Landau damping parameter (see Tables [T] 
& E|) suggests that very little of the turbulent en- 
ergy is dissipated via collisionless damping at large 
scales. In fact, even in the limit of zero conductiv- 
ity (k-L — * oo ) and thus no collisionless damping 
(the double-adiabatic limit) we find roughly the 
same energetics in our shearing box simulations. 
These results imply that heating due to work done 
by anisotropic stress dominates collisionless damp- 
ing on large scales in RIAFs. Collisionless damp- 
ing may become important on smaller scales (that 
we cannot resolve) where the turbulence can be 
better approximated as a superposition of linear 
modes (on large scales the fluctuating magnetic 
field is much larger than the mean field). Indeed, 
with higher resolution it may turn out that instead 
of being lost to grid-scale averaging, most of the 
kinetic and magnetic energy is damped via colli- 



7 Heating by Alfvenic turbulence is not expected to be impor- 
tant until very small scales of order the ion Larmor radius, 
scales that are not resolved by our MRI simulations. The 
fast mode is also strongly damped by collisionless damping, 
but is not likely to be as efficiently excited by the weakly 
compressible MRI turbulence. 



sionless damping at intermediate scales. In addi- 
tion, more accurate fully kinetic treatments may 
be needed to quantify the heating by collisionless 
damping on large scales in MRI simulations. 

Our results on viscous heating imply radiative 
efficiencies of > 0.5% for M > lQ- 4 M Ed d (Fig. 
|6|) . With such significant radiative efficiencies, the 
low luminosities of many accreting black holes can 
only be understood if most of the mass supplied 
to the flow at large radii does not reach the cen- 
tral object. Our results are thus consistent with 
global numerical simulations of magnetized RI- 
AFs, which show that only a small fraction of mass 
is actually accreted (e.g 



2001 



( Stone fc Pringle 
200lh . 

To give a concrete example of the implica- 
tions of our results, the Bondi accretion rate for 
Sgr A* in the Galactic Center is est imated to 
be ~ 3 x 10" 6 - lO^Mgyr" 1 (e.g., IQuataert 
120041: ICuadra et al. Il2006h . which corresponds to 
M/M^dd ~ 10~ 4 . Our results predict a radiative 
efficiency of « 0.005 - 0.03 for this M, with the 
uncertainty in the radiative efficiency arising from 
uncertainties in the exact rate of electron viscous 
heating. If gas were to accrete at the Bondi rate, 
a radiative efficiency of 0.01 implies a luminosity 
~ 10 3 times larger than what is observed from 
Sgr A*, implying that the accretion rate must be 
well below the Bondi rate. Requiring L ea 10 36 
ergs s _1 as is observed, our predicted radiative 
efficiencies (Fig. [5]) imply an accretion rate of 
M ~ 10~ 7 — lO _8 M0yr _1 where the range ac- 
counts for the range of radiative efficiencies in Fig- 
ure[6] Measurements of the rotation measure from 
Sgr A* are consistent with such a low accretion 

VLBI 
15R S 



rate ([Marrone et al. 1120071 ). In addition 



interferometr y of Sgr A* finds a size of 



at 100 GHz (|Bower et al. 112004 120061 ). The cor- 
responding brightness temperature is rs 3 x 10 10 
K. At low M, our models for the electron temper- 
ature at w 15i?s ar e in reasonable agreement with 
this result (Fig. |SJ). 
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Table 1 

Heating Diagnostics for Simulations with initially Cold Electrons (T*/T e = 10 4 at 5 

orbits) 





N x xN y x N z 




a M h 


CtR 


<#i/fipo d 






qjf V /np s 


0.5/fc 


27 x 59 x 27 


0.12 


0.097 


0.048 


0.18 


0.049 


0.14 


0.074 


0.5/Sz 


54 x 118 x 54 


0.21 


0.26 


0.083 


0.31 


0.06 


0.38 


0.13 


0.25/5z 


54 x 118 x 54 


0.18 


0.21 


0.073 


0.27 


0.043 


0.3 


0.11 


0.125/ 6z 


27 x 59 x 27 


0.16 


0.21 


0.084 


0.25 


0.04 


0.32 


0.13 


0.125/ 6z 


54 x 118 x 54 


0.15 


0.20 


0.060 


0.23 


0.049 


0.31 


0.094 


2/5z 


27 x 59 x 27 


0.16 


0.15 


0.073 


0.23 


0.016 


0.22 


0.11 


2/5z 


54 x 118 x 54 


0.19 


0.22 


0.071 


0.28 


0.039 


0.30 


0.11 


oo 


27 x 59 x 27 


0.19 


0.20 


0.086 


0.27 


0.013 


0.25 


0.13 


00 


54 x 118 x 54 


0.18 


0.22 


0.067 


0.27 


0.026 


0.27 


0.1 



a a J 4 = ({(p\\ — p±)(B x B y /B 2 )/po)), (()) denotes time and volume average in the turbulent state 
(from 5 to 20 orbits); po is the initial pressure. 

h a M - {{-B X B V /A%p )) 

c a R = {{pV x SV y /po)), where SV y = V y + (3/2)Qx 
d Qvi = 1Vi,i + qvi.e, where q vi , s is given in eq. [6] 
c Qv2 = 1V2,i + qv2, e , where q V 2,s is given in eq. [7] 
f Numerical loss of magnetic energy at the grid-scale 
s Numerical loss of kinetic energy at the grid-scale 



Table 2 

Heating Diagnostics for Different Ti/T e 



T t /T e a 


k L 


OLA 




OCR 


4 M /^o b 


Qv2,i/ n Po h 








Qnv/^Po 


1000 


0.5/Sz 


0.16 


0.17 


0.077 


0.16 


0.0037 


5.9 x 10~ 4 


1.6 X 10~ 4 


0.26 


0.12 


1000 


oo 


0.15 


0.13 


0.065 


0.16 


-0.0014 


8.2 x 10~ 4 


2.9 x 10~ 4 


0.17 


0.1 


100 


0.5/Sz 


0.13 


0.11 


0.057 


0.16 


0.0039 


0.003 


3.1 x 10~ 4 


0.17 


0.088 


100 


00 


0.18 


0.17 


0.079 


0.16 


-0.0024 


0.0042 


6.9 x 10~ 4 


0.23 


0.12 


10 


0.5/Sz 


0.16 


0.15 


0.071 


0.16 


0.0027 


0.01 


3.7 x 10~ 4 


0.23 


0.11 


10 


00 


0.15 


0.12 


0.061 


0.16 


-0.0055 


0.014 


0.0012 


0.16 


0.09 


1 


0.5/Sz 


0.17 


0.13 


0.07 


0.16 


-0.0015 


0.03 


2.1 x 10~ 4 


0.2 


0.11 


1 


oo 


0.18 


0.13 


0.068 


0.16 


-0.018 


0.035 


0.0056 


0.17 


0.1 



a The simulation is restarted after 5 orbits with this initial temperature ratio. N x x N y x N z = 27 X 59 x 27. All 
quantities have the same definitions as in Table [TJ 

b Because the electron and ion temperatures change significantly from 5 to 20 orbits, the heating rates are only 
averaged from 5 to 5.1 orbits so that the temperature is roughly equal to the value it is initialized at (column 1). 
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Fig. 1. — The increase in total energy normalized 
to the initial pressure po (solid line) in shearing 
box simulations of the MRI, compared to the work 
done by the stresses at the boundaries (dashed 
line). For a conservative calculation these two 
curves would agree identically. This is not the 
case because ZEUS is non-conservative. Nonethe- 
less, our method of improving energy co nservation 



in ZEUS (based on iTurner et al. 1 120031 ) conserves 



energy to better than ~ 5% over 15 orbits. 
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Fig. 3. — Volume averaged pressure anisotropy 
(AirAp/ B 2 ; Ap = p|| — p±) for electrons (thin solid 
line) and ions (thick solid line), compared to the 
anisotropy thresholds for the ion cyclotron (thick 
dashed line), mirror (thick dot-dashed line), and 
electron whistler (thin dashed line) instabilities. 
The ion to electron temperature ratio is set to 10 
after the initial saturation of the MRI at 5 orbits. 



Fig. 2. — Different contributions to plasma heat- 
ing in shearing box simulations of the MRI for 
initially cold electrons (Tj/T e = 10 4 at 5 orbits). 
Left: resolution of 27 x 59 x 27 resolution with 
Ul = 0-5/Sz, Right: resolution of 54x1 18x54 with 
kh = 0.25/(5z. qt is the heating due to anisotropic 
viscosity at large scales (electron + ion, but pri- 
marily ion because 7$ 3> T e ); q~^ R , the numerical 
loss of magnetic energy; q^ v , the numerical loss 
of kinetic energy. Viscous heating (qy) accounts 
for ~ 50% of the heating in our simulations. 
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Fig. 4. — The ratio of the volume averaged ion 
and electron temperatures as a function of time 
for models where the ion and electron tempera- 
ture ratios are set to 10 4 , 10 3 , 10 2 , 10, and 1 at 
the end of 5 orbits (top to bottom) . Solid lines are 
for kh = 0.5/Sz and dashed lines are for ki, = oo 
(the CGL limit). The electron heating is some- 
what larger in the CGL limit (see also Fig. [5])- 
Also shown (dot dashed line) is a calculation with 
Ti/T e = 10 3 at 5 orbits for k L = 0.5/Sz in which 
viscous heating of electrons is artificially turned 
off by setting qy e = 0. The electron temperature 
remains nearly constant but the ions are heated 
by viscosity, resulting in an increase in /T e with 
time. This calculation demonstrates that our sim- 
ulations can readily sustain large temperature dif- 
ferences between electrons and ions. 




Fig. 5. — The ion to electron heating ratio as a 
function of the temperature ratio for: k^ = 0.5/Sz, 
resolution 27 x 59 x 27 (+); k L = 0.25/£z, res- 
olution 54 x 118 x 54 (A); and the CGL limit, 
resolution 27 x 59 x 27 (o). The ion to electron 
heating ratio is slightly larger (less than a fac- 
tor of 2) with thermal conduction. The results 
with double the resolution are comparable to the 
standard resolution showing that the results are 
reasonably converged. The three solid lines corre- 
spond to q+/q+ = 6 (T,/^) 1 / 2 , 3 (T l /T e ) 1 /2 j and 
1.5 (Ti/Te) 1 / 2 (from top to bottom). These esti- 
mates of the electron heating rate, based on our 
analytic prediction from §2.2, approximate the nu- 
merical simulations reasonably well; they are used 
to calculate the radiative efficiency in §4 and Fig- 
ure [6] For this plot, the heating and temperature 
ratios are averaged over 0.1 orbits in a number of 
simulations (so that Ti/T e is roughly constant over 
the averaging). 
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Fig. 6. — Radiative efficiency as a function of 
mass accretion rate in one dimensional models 
of RIAFs (based on the cal culations described in 



Quataert &; Naravan 1999f ). These models 



sume a = 0.1 and /? = 10 in determining the 
density and magnetic field strength in the accre- 
tion flow. Since most of the radiation is produced 
close to the black hole, M can be interpreted as 
the accretion rate in the inner ~ 10i?s. The solid 
line corresponds to S = q+/qf = 0.33(T e /T l ) 1 / 2 , 
and the dashed lines to S = 0.66(T e /T l ) 1/2 (up- 
per) and S = 0.17(T e /Ti) 1/2 (lower). At high M, 
Coulomb collisions transfer energy from the ions 
to the electrons and so the radiative efficiency be- 
comes relatively independent of the assumptions 
regarding 5. At very low M, the electron cooling 
time is longer than the inflow time and the radia- 
tive efficiency decreases in all models. 
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Fig. 7. — Electron heating rate (volume averaged 
from 5 to 20 orbits) due to anisotropic viscosity as 
a function of the electron collisionality for simula- 
tions with no heat conduction (k^ — ► oo). The ion 
to electron temperature ratio is initialized to 10 3 
(*), 10 2 (+), 10 (V), & 1 (o) after 5 orbits. Viscous 
heating of electrons is suppressed at high electron 
collisionality (high v e ); the critical value of v e /£l in 
the numerical calculations is in reasonable agree- 
ment with the analytic estimate in equation ()15[) . 
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Fig. 8. — Electron temperature as a function of 
radius for our onc-dimcnsional RIAF models with 
a = 0.1, /3 = 10, and 5 = 0.33(T e /Tj) 1/2 . The 
solid lines correspond to log[M/MEdd] = -2, -2.5, 
-3, -3.5, -4, -4.5, -5, -5.5, -6 (from bottom to top). 
The electron temperature is independent of M for 
even lower accretion rates. The dotted line shows 
the electron temperature for S = 0.66(T e /Ti) 1 / 2 
and log[M/MEdd] = — 6. To demonstrate that 
at high M the electron temperature is relatively 
independent of our model assumptions, the dashed 
line shows T e (r) for M = 10~ 2 M E dd and 5 = 10~ 3 , 
while the dot-dashed line is for M = 10 _2 MEdd 
and S = OMiTe/Ti) 1 / 2 . 



